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Abstract 

The associated production of a VP boson with a jet originating from either a light 
parton or heavy-flavor quark is studied in the forward region using proton-proton 
collisions. The analysis uses data corresponding to integrated luminosities of 1.0 
and 2.0 fb~^ collected with the LHCb detector at center-of-mass energies of 7 and 
8 TeV, respectively. The VP bosons are reconstructed using the VP —>■ decay and 
muons with a transverse momentum, p^', larger than 20GeV in the pseudorapidity 
range 2.0 < p < 4.5. The partons are reconstructed as jets with px > 20GeV and 
2.2 < p < 4.2. The sum of the muon and jet momenta must satisfy px > 20 GeV. The 
fraction of VP-|-jet events that originate from beauty and charm quarks is measured, 
along with the charge asymmetries of the VP-l-6 and VP-|-c production cross sections. 
The ratio of the VP-Vjet to Z-|-jet production cross sections is also measured using 
the Z —7- pp decay. All results are in agreement with Standard Model predictions. 


Phys. Rev. D. 92 (2015) 052001 

© CERN on behalf of the LHCb collaboration, license CC-BY-4.0 


f Authors are listed at the end of this article. 




11 



1 Introduction 


Measurements of W^+jet production in hadron collisions provide important tests of the 
Standard Model (SM), especially of perturbative quantum chromodynamics (QCD) in the 
presence of heavy-flavor quarks. Such measurements are also sensitive probes of the parton 
distribution functions (PDFs) of the proton. The ratio of the IT+jet to Z-|-jet production 
cross sections is a test of perturbative QCD methods and constrains the light-parton PDFs 
of the proton. 

The jet produced in association with the W boson may originate either from a b quark 
(W + b), c quark (W + c) or light parton. Several processes contribute to the W+ b and 
W + c hnal states at next-to-leading order (NLO) in perturbative QCD. The dominant 
mechanism for W+c production is gs ^Wc, but there are also important contributions 
from gs —)■ Wcg, gg —)■ Wes, and qq ^ Wee [^. Therefore, measuring the ratio of the 
W+ c to IT-1-jet production cross sections in the forward region at the LHC provides 
important constraints on the s quark PDF at momentum transfers of Q ~ 100 GeV 
(c = 1 throughout this article) and momentum fractions down to a; ~ 10“®. Previous 
measurements of the proton s quark PDF were primarily based on deep inelastic scattering 
experiments with Q ~ 1 GeV and x values (P(O.l) JH- The W+c cross section has been 
measured at the Tevatron and at the LHG ilM in the central region. 

In the so-called four-flavor scheme, theoretical calculations are performed considering 
only the four lightest quarks in the proton 11 . Production of W+b proceeds via qq ^ Wg 
with g ^ bh a,t leading order. If the b quark content of the proton is considered, i.e. 

The 
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the hve-flavor scheme, then single-6 production via qb —)■ Wbq also contributes 
ratio of the lT-l-6 to lT-|-jet cross sections thus places constraints both on the intrinsic 6 
quark content of the proton and the probability of gluon splitting into 66 pairs. The lT-l-6 
cross section has been measured in the central region at the Tevatron 
LHG 
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13 14 and at the 


LHCb has measured the cross sections for inclusive W and Z production in proton¬ 
providing precision 


proton (pp) collisions at center-of-mass energy ^/s = 7TeV 16 19 


tests of the SM in the forward region. Additionally, measurements of the Z-|-jet and Z+b 
cross sections have been made 20,21 . In this article, the associated production of a bF 


boson with a jet originating from either a light parton or a heavy-flavor quark is studied 
using pp collisions at center-of-mass energies of 7 and 8TeV. The production of the lT-l-6 
hnal state via top quark decay is not included in the signal dehnition in this analysis, but 
is reported separately in Ref. 
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A comprehensive approach is taken, where the inclusive H^-|-jet, W + b and W + c 
contributions are measured simultaneously, rather than split across multiple measurements 
as in Refs. [l0| [I^ [23}j^ . The identiheation of c jets, in conjunction with 6 jets, 
is performed using the tagging algorithm described in Ref. |27 , which improves upon 


previous c-tagging methods where muons or exclusive decays were required to identify the 
jet [^[^. For each center-of-mass energy, the following production cross section ratios are 
measured: a(Wb)/a(Wj), a{Wc) j( t{W g), (j{W^g) j(j{Zg), (j{W~g) ja{Zg), A{Wb), and 
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A{Wc), where 


( 1 ) 


AiWq) 


a{W+q) - a{W-q) 
(j{W+q) + a{W-qy 


The analysis is performed using the hh —)■ /xz/ decay and jets clustered with the anti-fcx 
algorithm using a distance parameter R = 0.5. The following hducial requirements 
are applied: both the muon and the jet must have momentum transverse to the beam, px, 
greater than 20 GeV; the pseudorapidity of the muon must fall within 2.0 < ? 7 (/i) < 4.5; the 
jet pseudorapidity must satisfy 2.2 < rj{j) < 4.2; the muon and jet must be separated by 
> 0-5, where AR = ■>/and Ari{A(j)) is the difference in pseudorapidity 
(azimuthal angle) between the muon and jet momenta; and the transverse component of 
the sum of the muon and jet momenta must satisfy Pt(/^ + j) = i.PifJ') + p(j))t > ^0 GeV. 
All results reported in this article are for within this hducial region, i.e. no extrapolation 
outside of this region is performed. 

The article is organized as follows: the detector, data sample and simulation are 
described in Sect. the event selection is given in Sect. the signal yields are determined 
in Sect.lil the systematic uncertainties are outlined in Sect. [51 and the results are presented 
in Sect. El 


2 The LHCb detector and data set 


The LHGb detector |29|30l is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < p < 5, designed for the study of particles containing b or c quarks. The 
detector includes a high-precision tracking system consisting of a silicon-strip vertex 
detector surrounding the pp interaction region j^, a large-area silicon-strip detector 
located upstream of a dipole magnet with a bending power of about 4Tm, and three 
stations of silicon-strip detectors and straw drift tubes 
magnet. 
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placed downstream of the 
The tracking system provides a measurement of momentum, p, of charged 


particles with a relative uncertainty that varies from 0.5% at low momentum to 1.0% at 
200 GeV. The minimum distance of a track to a primary vertex, the impact parameter, 
is measured with a resolution of (15 -|- 29 /pt) M-ni; with px in GeV. Different types of 
charged hadrons are distinguished using information from two ring-imaging Cherenkov 
detectors. Photons, electrons and hadrons are identified by a calorimeter system consisting 
of scintillating-pad and preshower detectors, an electromagnetic calorimeter and a hadronic 
calorimeter. The electromagnetic and hadronic calorimeters have energy resolutions of 
a{E)/E = 10%/VE © 1% and a{E)/E = 69%/'/E © 9% (with E in GeV), respectively. 
Muons are identihed by a system composed of alternating layers of iron and multiwire 
proportional chambers 
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The trigger 34 consists of a hardware stage, based on information from the calorimeter 


and muon systems, followed by a software stage, which applies a full event reconstruction. 
This analysis requires at least one muon candidate that satishes the trigger requirement 
of px > 10 GeV. Global event cuts (GECs), which prevent high-occupancy events from 
dominating the processing time of the software trigger, are also applied and have an 
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efficiency of about 90% for W+]et and Z+jet events. 

Two sets of pp collision data collected with the LHCb detector are used: data collected 
during 2011 at s/s = 7TeV, corresponding to an integrated luminosity of 1.0 fb“^, and 
data collected during 2012 at a/s = 8 TeV, corresponding to an integrated luminosity of 
2.0 fb“/ Simulated pp collisions, used to study the detector response, to define the event 
selection and to validate data-driven techniques, are generated using Pythia 35 with 
an LHCb configuration 
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Decays of hadronic particles are described by EvtGen 

The interaction of 


in which final-state radiation (FSR) is generated using Photos 38 


the generated particles with the detector, and its response, are implemented using the 
Geant4 toolkit [Ml as described in Ref. (40 . 


Results are compared with theoretical calculations at NLO using MCFM and 
the CTIO PDF set |^. The theoretical uncertainty is a combination of PDF, scale, 
and strong-coupling (a^) uncertainties. The PDF and scale uncertainties are evaluated 
following Refs. and [^, respectively. The as uncertainty is evaluated as the envelope 
obtained using as{Mz) G [0.117,0.118,0.119] in the theory calculations. 


3 Event selection 


The signature for bF-l-jet events is an isolated high-px muon and a well-separated jet, both 
produced in the same pp interaction. Muon candidates are identified with tracks that have 
associated hits in the muon system. The muon candidate must have Pt(f) > 20 GeV and 
pseudorapidity within 2.0 < p(/i) < 4.5. Background muons from kF —)■ r —)■ p decays or 
semileptonic decays of heavy-flavor hadrons are suppressed by requiring the muon impact 
parameter to be less than 0.04 mm |16 . Background from high-momentum kaons and 


pions that enter the muon system and are misidentified as muons, is reduced by requiring 
that the sum of the energy of the associated electromagnetic and hadronic calorimeter 
deposits does not exceed 4% of the momentum of the muon candidate. 

Jets are clustered using the anti-Zcx algorithm with a distance parameter R = 0.5, as 
implemented in Fastjet 44 . Information from all the detector subsystems is used to 


create charged and neutral particle inputs to the jet-clustering algorithm using a particle 
flow approach |^. During 2011 and 2012, LHCb collected data with a mean number 
of pp collisions per beam crossing of about 1.7. To reduce contamination from multiple 
pp interactions, charged particles reconstructed within the vertex detector may only be 
clustered into a jet if they are associated with the same pp collision. 

Signal events are selected by requiring a muon candidate and at least one jet with 
AR{p,j) > 0.5. For each event the highest-px muon candidate that satisfies the trigger 
requirements is selected, along with the highest-px jet from the same pp collision. The 
high-px muon candidate is not removed from the anti-fcx inputs and so is clustered into 
a jet. This jet, referred to as the muon jet and denoted as jji, is used to discriminate 
between kF-|-jet and dijet events. The requirement PT{j^l+j) > 20 GeV is made to suppress 
dijet backgrounds, which are well balanced in px, unlike kF-|-jet events where there is 
undetected energy from the neutrino. Furthermore, the distribution of the fractional muon 
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candidate px within the muon jet, /VtU^ l)i is used to separate vector bosons from 

jets. For vector-boson production, this ratio deviates from unity only due to muon FSR, 
activity from the underlying event, or from neutral-particle production in a separate pp 
collision, whereas for jet production this ratio is driven to smaller values by the presence 
of additional radiation produced in association with the muon candidate. 

Events with a second, oppositely charged, muon candidate from the same pp collision 
are vetoed. However, when the dimuon invariant mass is in the range 60 < M{pi^p~) < 
120 GeV, such events are selected as Z-|-jet candidates and the PT^{i^i+i) requirement is not 
applied. Two Z-|-jet data samples are selected at each center-of-mass energy: a data sample 
where only the /i^ is required to satisfy the trigger requirements and one where only the 
/i“ is required to satisfy them. The first sample is used to measure a{W^i)/a{Zi), while 
the second is used for a{W~j)/a{Zi). This strategy leads to approximate cancellation of 
the uncertainty in the trigger efficiency in the measurement of these ratios. 

The reconstructed jets must have PtU) > 20 GeV and 2.2 < p(j) < 4.2. The reduced 
p{j) acceptance ensures nearly uniform jet reconstruction and heavy-flavor tagging efficien¬ 
cies. The momentum of a reconstructed jet is scaled to obtain an unbiased estimate of the 
true jet momentum. The scaling factor, typically between 0.9 and 1.1, is determined from 
simulation and depends on the jet px and p, the fraction of the jet transverse momentum 
measured with the tracking systems, and the number of pp interactions in the event. No 
scaling is applied to the momentum of the muon jet. Migration of events in and out of 
the jet px hducial region due to the detector response is corrected for by an unfolding 
technique. Data-driven methods are used to obtain the unfolding matrix, with the resulting 
corrections to the measurements presented in this article being at the percent level. 

The jets are identified, or tagged, as originating from the hadronization of a heavy-flavor 
quark by the presence of a secondary vertex (SV) with Ai? < 0.5 between the jet axis and 
the SV direction of flight, dehned by the vector from the pp interaction point to the SV 
position. Two boosted decision trees (BDTs) |4^[4^ , 'SDT{hc\udsg) and BDT(&|c), trained 
on the characteristics of the SV and the jet, are used to separate heavy-flavor jets from 
light-parton jets, and to separate b jets from c jets. The two-dimensional distribution of 
the BDT response observed in data is fitted to obtain the SV-tagged b, c and light-parton 
jet yields. The SV-tagger algorithm is detailed in Ref. 27 , where the heavy-flavor tagging 


efficiencies and light-parton mistag probabilities are measured in data. 


4 Background determination 

Contributions from six processes are considered in the fF-|-jet data sample: bF-|-jet signal 
events; Z-|-jet events where one muon is not reconstructed; top quark events producing a 
bF-|-jet hnal state; Z —rr events where one r lepton decays to a muon and the other 
decays hadronically; QCD dijet events; and vector boson pair production. Simulations 
based on NLO predictions show that the last contribution is negligible. 

The signal yields are obtained for each muon charge and center-of-mass energy in¬ 
dependently. The Pt{p)!PT:{ i distribution is htted to determine the bF-|-jet yield of 
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each data sample. To determine the W + h and W + c yields, the snbset of candidates 
with an SV-tagged jet is binned according to /VtH^ l)■ each PT(h)/PT(j/i) bin, the 

two-dimensional SV-tagger BDT-response distributions are htted to determine the yields 
of 6-tagged and c-tagged jets, which are used to form the V^ip)!V tUij) distributions for 
candidates with 6-tagged and c-tagged jets. These ^^) distributions are htted 

to determine the SV-tagged W + h and W + c yields. Finally, to obtain a{Wh)/a{Wj) 
and a{Wc)/a{Wi), the jet-tagging efficiencies of etag(&) ~ 65% and etag(c) ~ 25% are 
accounted for. In all hts performed in this analysis, the templates are histograms with 
hxed shapes. 

The PT(/^)/PT(jji) distributions are shown in Fig. (in this and subsequent hgures the 
pull represents the difference between the data and the ht, in units of standard deviations). 
The W boson yields are determined by performing binned extended-maximum-likelihood 
hts to these distributions with the following components: 

• The W boson template is obtained by correcting the !V'tU fi) distribution 
observed in Z+jet events for small differences between W and Z decays derived from 
simulation. 


• The template for Z boson events where one muon is not reconstructed is obtained 
by correcting, using simulation, the ^^) distribution observed in fully 

reconstructed Z-|-jet events for small differences expected in partially reconstructed 
Z-|-jet events. The yield is hxed from the fully reconstructed Z-|-jet data sample, 
where simulation is used to obtain the probability that the muon is missed, either 
because it is out of acceptance or it is not reconstructed. 


• The templates for 6, c and light-parton jets are obtained using dijet-enriched data 
samples. These samples require pT(jVi+j) < lOGeV and, for the heavy-havor samples, 
either a stringent 6-tag or c-tag requirement on the associated jet. The templates 
are corrected for diherences in the Pt(%) spectra between the dijet-enriched and 
signal regions. The contributions of b, c and light-parton jets are each free to vary 
in the PT(h)/PT(%) 


The PT{p)/PT{jfi) hts determine the bF+jet yields, which include contributions from top 
quark and Z ^ tt production. The top quark and Z ^ tt contributions cannot be 
separated from VF+jet since their Pt{p)/PtU fi) distributions are nearly identical to that 
of VF+jet events. The subtraction of these backgrounds is described below. 

The yields of events with IF bosons associated with 6-tagged and c-tagged jets are 
obtained by htting the two-dimensional SV-tagger BDT-response distributions for ^/s = 7 
and 8TeV and for each muon charge separately in bins of P'lip)/pT{i^l)■ The SV-tagger 
BDT templates used in this analysis are obtained from the data samples enriched in b and 
c jets used in Ref. 27 . As a consistency check, the two-dimensional BDT distributions 


are htted using templates from simulation; the yields shift only by a few percent. Figure 
shows the BDT distributions combining all data in the most sensitive region, IF-|-jet events 
with Pt(f)/Pt(%) > 0.9. This is the region where the muon carries a large fraction of the 
muon-jet momentum and is, therefore, highly isolated. Figure shows the distributions in 
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Figure 1: Distributions of with fits overlaid from (top) ^/s = 7TeV and (bottom) 

8TeV data for (left) and (right) p~. 


a dijet dominated region (0.5 < /PtU ii) < 0.6). In the dijet region the majority of 

SV-tagged jets associated with the high-px muon candidate are found to be b jets. This is 
due to the large semileptonic branching fraction of b hadrons. In the IT+jet signal region 
there are significant contributions from both b and c jets. 

As a consistency check, the 6, c, and light-parton yields are obtained in the 
/PtU fi) > 0.9 signal region from a ht nsing only two of the BDT inpnts, both 
of which rely only on basic SV properties, the track multiplicity and the corrected mass, 
which is defined as 

Mcor = \l + |p|2gin2 Q _|_ 1^ sin 51^ ( 2 ) 

where M and p are the invariant mass and momentum of the particles that form the SV, and 
9 is the angle between p and the flight direction. The corrected mass, which is the minimum 
mass for a long-lived hadron whose trajectory is consistent with the flight direction, peaks 
near the D meson mass for c jets and conseqnently provides excellent discrimination 
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Figure 2: (top left) Two-dimensional SV-tag BDT distribution and (top right) fit for events in 
the subsample with PT{^J')/PT{j^l) > 0.9, projected onto the (bottom left) BDT{bc\udsg) and 
(bottom right) BDT(6|c) axes. Combined data for ^/s = 7 and 8 TeV for both muon charges are 
shown. 


against other jet types. The SV track multiplicity identihes b jets well, since 6-hadron 
decays typically produce many displaced tracks. In Fig. the distributions of Mcor and 
SV track multiplicity for a subsample of SV-tagged events with BDT{bc\udsg) > 0.2 (see 
Fig. 1^ are htted simultaneously. The templates used in these hts are obtained from data 
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Figure 3: (top left) Two-dimensional SV-tag BDT distribution and (top right) fit for events in 
the subsample with 0.5 < PTib^)/P'rii^i) < 0.6, projected onto the (bottom left) BDT( 6 c|u(is 5 f) 
and (bottom right) BDT( 6 |c) axes. Combined data for ^/s = 7 and 8 TeV for both muon charges 
are shown. 


in the same manner as the SV-tagger BDT templates. After correcting for the efficiency 
of requiring BDT{bc\udsg) > 0.2, the b and c yields determined from the hts to Mcor and 
SV track multiplicity and from the two-dimensional BDT hts are consistent. The mistag 
probability for hT-|-hght-parton events in this sample is found to be approximately 0.3%, 
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Figure 4: Projections of simultaneous fits of (left) Mcor and (right) SV track multiplicity for the 
SV-tagged subsample with BDT{bc\udsg) > 0.2 and > 0.9. The highest Mcor bin 

includes candidates with Mcor > 10 GeV. Combined data for y/s = 7 and 8 TeV for both muon 
charges are shown. 


which agrees with the value obtained from simulation. 

From the SV-tagger BDT hts, the b and c yields are obtained in bins of y/s, muon 
charge, and PT{p)/PT{jiJ,)- The Pt:{p)/ptU^) distribntions for mnons associated with 
6 -tagged and c-tagged jets are shown in Figs. [^and|^ These distributions are htted to 
determine the hF+6 and W-hc hnal-state yields as in the inclnsive hF+jet sample. The 
Z-hb and Z+c yields are obtained by fitting the SV-tagger BDT distributions in the fully 
reconstructed Z-|-jet data samples and then correcting for the missed-muon probability. 
The fits are shown in Figs. and for each muon charge and center-of-mass energy. The 
yields obtained still include contributions from top quark production and Z —)■ rr. 

The Z ^ TT background, where one r lepton decays into a muon and the other into a 
hadronic jet, contaminates the IT-fc sample dne to the similarity of the c-hadron and r 
lepton masses. The Pt(SV)/pt(j) distribution, where pt(SV) is the transverse momentum 
of the particles that form the SV, is used to discriminate between c and r jets, since SVs 
prodnced from r decays nsnally carry a larger fraction of the jet energy than SVs from 
c-hadron decays. Figure shows fits to the PT(SV)/pT(j) distributions observed in data 
where the b and light-parton yields are fixed using the results of BDT fits performed on the 
data samples. A requirement of BDT{bc\udsg) > 0.2 is applied to this sample to remove 
the majority of SV-tagged light-parton jets while retaining 90% of b, c and r jets. The 
only free parameter in these fits is the fraction of jets identified as charm in the SV-tagger 
BDT hts that originate from r leptons. The Pt(SV)/pt(j) templates are obtained from 
simnlation. The Z —)■ rr yields are consistent with SM expectations and are abont 25 
times smaller than the W+c yields. These results are extrapolated to the inclusive sample 
using simnlation. 

The top quark background is determined in the dedicated analysis of Ref. 22 , where 
a rednced hducial region is nsed to enrich the relative top quark content. The yields 
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Figure 5: Fits to /PT{jp) distributions for 6-tagged data samples for y/s = 7 and 8 TeV. 


and charge asymmetries of the hF + 6 hnal state as functions of + b) are used to 
discriminate between W+h and top quark production. The results obtained in Ref. 
are consistent with SM expectations and are extrapolated to the hducial region of this 
analysis using simulation based on NLO calculations. The extrapolated top quark yields 
are subtracted from the observed number oiW+h candidates to obtain the signal yields. 
Top quark production is found to be responsible for about 1/3 of events that contain a W 
boson and 6 jet. A summary of all signal yields is given in Table 


5 Systematic uncertainties 

A summary of the relative systematic uncertainties separated by source for each mea¬ 
surement is provided in Table A detailed description of each contribution is given 
below. 

The pt distributions of muons from W and Z bosons produced in association with 
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Figure 6: Fits to PT(i^)/PT(j^) distributions for c-tagged data samples for y/s = 7 and 8 TeV. 


b, c and light-parton jets are nearly identical. This results in a negligible uncertainty 
from muon trigger and reconstruction efficiency on cross section ratios involving only W 
bosons. In the ratios a{W~^j)/a{Zj) and a{W~j)/a{Zj), the muon from the Z boson 
decay with the same charge as that from the W decay is required to satisfy the same 
trigger and selection requirements as the W boson muon, giving negligible uncertainty 
from the trigger and selection efficiency. The efficiency for reconstructing and selecting 
the additional muon from the Z boson decay is obtained from the data-driven studies of 
Ref. 17 . A further data-driven correction is applied to account for the higher occupancy 
in events with jets [^; a 2% systematic uncertainty is assigned to this correction. 

The GEC efficiency is obtained following Ref. 1^: an alternative dimuon trigger 


requirement with a looser GEC is used to determine the fraction of events that are rejected. 
The GEC efficiencies for all hnal states are found to be consistent within a statistical 
precision of 1%, which is assigned as a systematic uncertainty. As a further check, the 
number of jets per event reconstructed in association with W or Z bosons is compared 
and found to be consistent. 
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Figure 7: Fits to the PT(SV)/pT(j) distributions in (left) 7TeV and (right) 8 TeV data for 
candidates with PT(^)/ 7 'T(j/i) > 0.9 and BDT( 6 c|M(is 5 ) > 0.2. 


Table 1: Summary of signal yields. The two Zj yields denote the charge of the muon on which 
the trigger requirement is made. The Zj yields given are the numbers of candidates observed, 
while the W boson yields are obtained from fits. The yield due to top quark production is 
subtracted in these results. 



7TeV 

8 TeV 

Mode 

p+ 

/i 

/i+ 

/i 


2364 

2357 

6680 

6633 

Wj 

27400 ±500 

17 500 ±400 

70 700 ± 1100 

44 800 ± 800 

Wb-ia.g 

160 ± 31 

51 ±27 

400 ± 43 

236 ± 45 

ITc-tag 

295 ± 36 

338 ± 31 

795 ± 56 

802 ± 55 


The jet reconstruction efficiencies for heavy-flavor and light-parton jets in simulation 
are found to be consistent within 2 %, which is assigned as a systematic uncertainty for 
flavor-dependencies in the jet-reconstruction efficiency. The jet pt detector response 
is studied with a data sample enriched in h jets using SV tagging. The PT(SV)/pT(j) 
distribution observed in data is compared to templates obtained from simulation in bins 
of jet pt- The resolution and scale in simulation for each jet px bin are varied to find the 
best description of the data and to construct a data-driven unfolding matrix. The results 
obtained using this unfolding matrix are consistent with those obtained using a matrix 
determined by studies of px balance in Z-l-jet events 20 , where no heavy-flavor tagging is 
applied. The unfolding corrections are at the percent level and their statistical precision is 
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assigned as the uncertainty. 

The heavy-flavor tagging efficiencies are measured from data in Ref. [^, where a 10% 
uncertainty is assigned for b and c jets. The cross-check fits of Sect. using the corrected 
mass and track multiplicity, remove information associated with jet quantities, such as 
Pt, from the yield determination and produce yields consistent at the 5% level. This is 
assigned as the uncertainty for the SV-tagged yield determination. 

The W boson template for the PT{b'^)/PT{jii) distribution is derived from data, as 
described in Sect. The fit is repeated using variations of this template, e.g. using a 
template taken directly from simulation and using separate templates for W~^ and W~, 
to assess a systematic uncertainty. The dijet templates are obtained from data in a 
dijet-enriched region. The residual, small W boson contamination is subtracted using two 
methods: the W boson yield expected in the dijet-enriched region is taken from simulation; 
and the Pt{p)/PtU^ i) distribution in the dijet-enriched region is fitted to a parametric 
function to estimate the W boson yield. The difference in the W boson yields obtained 
using these two sets of dijet templates is at most 2%. The uncertainty on W/Z ratios due 
to the W boson and dijet templates is 4%. The uncertainty due to the W boson template 
cancels to good approximation in the measurements of a{Wh)/a{Wj) and a{Wc)/a{Wj)] 
however, the uncertainty due to the dijet templates is larger due to the enhanced dijet 
background levels. Variations of the dijet templates are considered, with 10% and 5% 
uncertainties assigned on a{Wh)/a{Wj) and a{Wc)/a{Wj). 

The systematic uncertainty from top quark production is taken from Ref. [^, while 
the systematic uncertainty from Z —)■ rr is evaluated by fitting the data using variations of 
the PT(SV)/pT(j) templates. All other electroweak backgrounds are found to be negligible 
from NLO predictions. All W ^ pu yields have a small contamination from W ^ t ^ p 
decays that cancels in all cross section ratios except for the W/Z ratios. A scaling factor 
of 0.975, obtained from simulation, is applied to the W boson yields. A 1% uncertainty is 
assigned to the scale factor, which is obtained from the difference between the correction 
factor from simulation and a data-driven study of this background for inclusive 
W ^ pu production. 

The trigger, reconstruction and selection requirements are consistent with being charge 
symmetric 16 , which results in negligible uncertainty on A{Wh) and AiWc). Unfolding of 
the jet pt detector response is performed independently for and W~ bosons, with the 
statistical uncertainties on the corrections to the charge asymmetries assigned as systematic 
uncertainties. The uncertainty on the WAh and WAc yields from the BDT templates 
is included in the charge asymmetry uncertainty due to the fact that the fractional jet 
content of the SV-tagged samples is charge dependent. The uncertainty on the charge 
asymmetries due to determination of the W boson yields is evaluated using an alternative 
method for obtaining the charge asymmetries. The raw charge asymmetry in the &-jet 
and c-jet yields in the Pt{p)/Pt:{ 3^l) > 0.9 region is obtained from the SV-tagger BDT 
fits. The Z-|-jet and dijet backgrounds are charge symmetric at the percent level and 
contribute at most to 20% of the events in this Pt{p) /P^Uii) region. Therefore, A{Wh) 
and A{Wc) are approximated by scaling the raw asymmetries by the inverse of the W 
boson purity in the Pt{p) /PT:{i^i) > 0.9 region. A small correction must also be applied 
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Table 2: Systematic uncertainties. Relative uncertainties are given for cross section ratios and 
absolute uncertainties for charge asymmetries. 


Source 

(7{Wb) 

a(Wj) 

a-{Wc) 

a{Wj) 

a{Wj) 

A(Wb) 

A(Wc) 

Muon trigger and selection 

— 

— 

2% 

— 

— 

GEC 

1 % 

1 % 

1% 

— 

— 

Jet reconstrnction 

2% 

2% 

— 

— 

— 

Jet pt 

2% 

2% 

1% 

0.02 

0.02 

(6, c)-tag efficiency 

10% 

10% 


— 

— 

SV-tag BDT templates 

5% 

5% 


0.02 

0.02 

PT{R)/PT{j^,) templates 

10% 

5% 

4% 

0.08 

0.03 

Top quark 

13% 

— 

— 

0.02 


Z ^ TT 

— 

3% 

— 

— 

— 

Other electroweak 

— 

— 

— 

— 

— 

W ^ T p, 

— 

— 

1 % 

— 

— 

Total 

20% 

13% 

5% 

0.09 

0.04 


Table 3: Summary of the results and SM predictions. For each measurement the first uncertainty 
is statistical, while the second is systematic. All results are reported within a fiducial region 
that requires a jet with px > 20GeV in the pseudorapidity range 2.2 < rj < 4.2, a muon with 
Pt > 20GeV in the pseudorapidity range 2.0 < ry < 4.5, PTifJ^ + j) > 20GeV, and AR{ij,,j) > 0.5. 
For Z+jet events both muons must fulfill the muon requirements and 60 < < 120 GeV; 

the Z+jet fiducial region does not require Pt(^ + i) > 20 GeV. 



Results 

7TeV 8TeV 

SM prediction 

7TeV 8TeV 

X 10^ 

cr(Wj) 

izi X10^ 

cr{Wj) 

0.66 ±0.13 ±0.13 

5.80 ±0.44 ±0.75 

0.78 ±0.08 ±0.16 

5.62 ±0.28 ±0.73 

n 74+0-1'' 

'+-0.13 

5 

+ -'J+-0.69 

n 77+0.18 

•O- ' '-0.13 
c qi +0.87 

A{Wb) 

aIwc) 

0.51 ±0.20 ±0.09 

-0.09 ±0.08 ±0.04 

0.27 ±0.13 ±0.09 
-0.01 ±0.05 ±0.04 

fi o'T+O.OS 
^ -0.03 
— 0 1 5+0-02 
"O--‘-"1-0.04 

n oq+O-03 
U*^°-0.03 

—n 14 + 0.02 

'0--‘-+-0.03 

a{W+j) 

<Zj) 

a(W-j) 

^{Zj) 

10.49 ±0.28 ±0.53 

6.61 ±0.19 ±0.33 

9.44 ±0.19 ±0.47 

6.02 ±0.13 ±0.30 

9 qn+0-28 

y.yu_Q 24 

5 70 + 0.21 
+ • '11-0.18 

q 48 + 0.16 

ir [r9+0.13 
+ •++-0.25 


to A{Wb) to account for top quark production. The difference between the asymmetries 
from this method and the nominal method is assigned as a systematic nncertainty from 
W boson signal determination. The nncertainty on A(Wb) due to top quark production 
is taken from Ref. 
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6 Results 


The results for y/s = 7 and 8 TeV are summarized in Table Each result is compared 
to SM predictions calculated at NLO using MCFM 


41 and the CTIO PDF set 42 


as 

described in Sect. Production of IT+jet events in the forward region requires a large 
imbalance in x of the initial partons. In the four-flavor scheme at leading order, W+ b 
production proceeds via qq —)■ Wg{hh), where the charge of the W boson has the same 
sign as that of the initial parton with larger x. Therefore, A{Wh) ~ -1-1/3 is predicted 
due to the valence quark content of the proton. The dominant mechanism for kP -t- c 
production is gs ^ Wc, which is charge symmetric assuming symmetric s and s quark 
PDFs. However, the Cabibbo-suppressed contribution from gd ^ Wc leads to a prediction 
of a small negative value for A(Wc). 

The a{Wb)/a{Wj) ratio in conjunction with the W+b charge asymmetry is consistent 
with MCFM calculations performed in the four-flavor scheme, where W+b production is 
primarily from gluon splitting. This scheme assumes no intrinsic b quark content in the 
proton. The data do not support a large contribution from intrinsic b quark content in 
the proton but the precision is not sufficient to rule out such a contribution at 0(10%). 
The ratio [a{Wb) + (T(top)]/cr(kFj) is measured to be 1.17 ± 0.13 (stat) ± 0.18 (syst)% at 
-y/s = 7 TeV and 1.29 ± 0.08 (stat) ± 0.19 (syst)% at y/s = 8 TeV, which agree with the 
NLO SM predictions of 1.23 ± 0.24% and 1.38 ± 0.26%, respectively. 

The a{Wc)/a{Wj) ratio is much larger than a{Wb)/a{Wj), which is consistent with 
Wc production from intrinsic s quark content of the proton. The measured charge 
asymmetry for lF-|-c is about 2a smaller than the predicted value obtained with CTIO, 
which assumes symmetric s and s quark PDFs. This could suggest a larger than expected 
contribution from scattering off of strange quarks or a charge asymmetry between s and 
s quarks in the proton. The ratio a{W~^j)/a{Zj) is consistent within la with NLO 
predictions, while the observed a{W~j)/a{Zj) ratio is higher than the predicted value by 
about 1.5a. 
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